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Abstract 

In  section  1  we  present  a  general  principle  for 
associating  nonlinear  equations  of  evolutions  with 
linear  operators  so  that  the  eigenvalues  of  the  linear 
operator  are  integrals  of  the  nonlinear  equation.   A 
striking  instance  of  such  a  procedure  is  the  discovery 
by  Gardner,  Miura  and  Kruskal  that  the  eigenvalues  of 
the  Schr^dinger  operator  are  integrals  of  the  Korteweg- 
de  Vries  equation. 

In  section  2  we  prove  the  simplest  case  of  a 
conjecture  of  Kruskal  and  Zabusky  concerning  the  existence 
of  double  wave  solutions  of  the  Korteweg-de  Vries  equation, 
i.e.  of  solutions  which  for   |t|   large  behave  as  the  super- 
position of  two  solitary  waves  travelling  at  different  speeds 
The  main  tool  used  is  the  first  of  a  remarkable  series   of 
integrals  discovered  by  Kruskal  and  Zabusky. 


Ill 


INTEGRALS  OF  NONLINEAR  EQUATIONS  OF 
EVOLUTION  AND  SOLITARY  WAVES 


In  this  paper  we  study  the  equation 

(1-1)  ^t  +  ^^x  +^xxx  =  0 

introduced  by  Korteweg  and  de  Vries  in  their  approximate  theory 
of  water  waves,   [3];  we  shall  refer  to  it  as  the  KdV  equation. 
Subsequently  the  KdV  equation  was  found  to  be  relevant  for  the 
description  of  hydromagnetic  waves  [ 2] ,   and  in  the  description 
of  acoustic  waves  in  an  anharmonic  crystal,  [8] .    (1.1)   is  a 
special  instance  of  a  nonlinear  evolution  equation  of  form 

(1.2)  u^  =  K(u) . 

We  shall  study   C   solutions  of  (1.1)  defined  for  all 
X  — ^  +  oo ,   together  with  all  their  x  derivatives.   It  is 
easy  to  show  that  such  solutions  are  uniquely  determined  by 
their  initial  values;  for  let   v  be  another  solution  of  (1.1): 

(1.1)  V,  +W+V    =0 

^      V  t        X      XXX 

Subtracting  this  from  (1.1)  and  denoting  u  -  v  by  w  we 
obtain  the  linear  equation 

^t  +  ^^x  +  ^x  +  ^xxx  =  ° 
for  w.   Multiplying  by  w  and  integrating  with  respect  to  x 


over   (-00,00)   we  obtain,  after  integration  by  parts,  and 
using  the  fact  that  w  and  w    tend  to  zero  as   x  — >  +  od 
the  relation 

(1.3)  ^  |/-'d-  +/^"-  "  l^x^^'^^  =  0  • 

Denoting   1/2  /w  dx   by  E(t)   and  Max ] 2v^-u^ |   by  m  we 
obtain  from  (1.3)  the  inequality 

A  E(t)  <  mE(t). 
This  differential  inequality  implies  that 

(1.4)  E(t)  <  £(0)6""^ 

which  implies  in  particular  that  if  E(0)   is  zero  then  so 
is   E(t)  -  and  thereby  w  -  for  all  t.   Furthermore, 
Sjoberg  has  shown,  [6],    that  (1.1)  has  a  solution  with 
arbitrarily  prescribed  initial  value   f(x)   provided  that 
f  is  smooth  enough  and  tends  to  zero  with  sufficient  rapidity 
as   |x|   tends  to  infinity. 

Equation  (1.1)  has  travelling  wave  solutions,  i.e.  of 
the  form  u(x,t)  =  s(x-ct),   c   being  the  speed  of  the  wave. 
To  see  this  substitute  into  (l.l): 

(1.5)  -  cs   +  ss   +  s    =  0  . 

^    -^  '  X       X      XXX 

Integrating  with  respect  to  x  and  imposing  the  boundary 
condition  that   s   and  its  derivatives  vanish  at  x  =  +  00 


we  get 

,2 
2~    '  "xx 


(1.6)         -  cs  +  -^  +  s   =  0, 


Multiplication  by   2s    and  one  more  Integration  gives 

(1.7)  -  cs^  +  J-  +  s^  =  0. 

From  this  relation   s   can  be  determined  explicitly: 

(1.8)  ■      s(x)  -  3c  sech^(^/^) 

Thus  we  see  that  (1.5)  has  a  solution  vanishing  at   x  =  +  oo 
for  every  positive  speed   c,   uniquely  determined  except  for 
a  shift  which  can  be  so  chosen  that  the  maximum  of   s   occurs 
at  X  =  0.   We  denote  this  normalized   s,   explicitly  described 
by  (1.8),  as_  s(x;c);   s(x,c)   is  symmetric  in  x,   decays 
exponentially  as   |x|  —^   oo ,   and   s(0,c)  =:  3c. 

On  account  of  its  shape   s   is  called  a  solitary  wave. 

For  linear  equations  it  often  happens  that  all  solutions 
can  be  obtained  as   superpositions  of  a  family  of  special 
solutions;  e.g.  all  solutions  of  linear  equations  with 
constant  coefficients  are  superpositions  of  exponential 
solutions.   For  nonlinear  equations  one  cannot  in  general 
form  new  solutions  out  of  old  and  so  special  families  of 
solutions  are  not  expected  to  play  any  special  role  in  the 
description  of  all  solutions.   It  was  therefore  very  surprising 
when  Kruskal  and  Zabusky  observed,  by  analyzing  numerically 


computed  solutions,  that  all  solutions  of  the  KdV  equation 
have  hidden  in  them  solitary  waves.  A  precise  formulation 
of  this  observation  is  as  follows: 

Let  u  be  any  solution  of  (l.l)  which  is  defined  for 
all   X  and   t   and  which  vanishes  at   x  =  +  ^  •   Then  there 
is  a  discrete  set  of  positive  numbers   c-,  ,  .  .  .,  C^^   -  called 
the  eigenspeeds  of  u  -   and  sets  of  phase  shifts   9^ 


such  that 


s(x-0i,c  .)   if  c  -  c  . 


(1.9)       lim  u(x+ct,t)  = 

^^  ±  °°  0     if  c  ;^  c  . 

J 

Note:  Consider  any  equation  of  evolution  (1.2)  which 
does  not  involve  x  and   t   explicitly,  i.e.  whose  set  of 
solutions  is  invariant  under  translation  with  respect  to 
X  and   t.   Suppose  u(x,t)   is  a  solution  of  such  an  equation, 
and  that  for  a  certain  value  of  c 

lim  u(x+ct,t) 

t— >  +  00 

exists  uniformly  on  compact  sets  in  x   space.   Clearly  the 
limit  is  a  travelling  wave  solution;  one  could  then  define 
the  eigenspeeds  of  a  solution  u  as  those  values  of   c   for 
which  the  above  limit  exists  and  is  different  from  zero.   It 
is  far  from  obvious  that  -  as  is  the  case  for  the  KdV  equation 
the  eigenspeeds  which  appear  in  the  limit   t  — >  oo   and 
t  — >  -  00   are  the  same. 


The  eigenspeeds   c  .   are  unequivocally  determined  by 
the  solution  u  under  consideration  and  thus  can  be  regarded 
as  functionals  of  solutions.   It  is  clear  from  the  definition 
that  if  u   is  translated  by  the  amounts   a  and  b   in  the 
X  and  t   direction,  the  eigenspeeds  remain  the  same  while 
the  phase  shifts  change  by  the  amount   a  -  cb. 

We  saw  earlier  that  solutions  are  uniquely  determined 

by  their  initial  values;  therefore  the  eigenspeeds  can  be  also 

regarded  as  functionals  of  the  initial  values.   It  follows 

from  translation  invariance  that  the  eigenspeeds  are. 

invariant  functionals  (also  called  Integrals);  that  is,  if 

f  and   f   denote  the  value  of  u  at  two  different  times, 

then   c.(f)  =  c.(f').    It  follows  similarly  that  the  difference 

© .  -  0~   is  an  integral. 
J     J 

This  analysis  shows  that  if  solutions  of  the  KdV 
equation  behave  for  large   t   as  indicated  in  (1.9)  then  the 
KdV  equation  has  an  infinity  of  integrals.   Indeed,  Kruskal, 
Gardner  and  Mlura,  [5],   succeeded  in  constructing  explicitly 
an  infinite  sequence  of  integrals;  the  simplest  of  these 
will  be  described  and  used  in  section  2.   We  present  now  a 
general  method  for  constructing  an  infinite  set  of  integrals 
for  equations  of  evolutions  (1.2). 

Let  (p   be  some  space  of  functions  chosen  so  that 
for  each  t,   u(t)   lies  in  (p.      Suppose  that  to  each  function 
u   in  C/p  we  can  associate  a  selfadjoint  operator  L  =  L 


over  some  Hilbert  space: 

(1.10)  ^  -^  ^u  V 

with  the  following  property:  That  if  u  changes  with   t 
subject  to  the  equation 

u^  =  K(u) 

the  operators   L{t),   which  also  change  with  t,      remain 
unitarily  equivalent.   If  this  is  the  case  then  the  eigenvalues 
of  L   constitute  a  set  of  integrals  for  the  equation  under 

U  D  T. 

consideration. 

The  unitary  equivalence  of  the  operators   L(t)   means 
that  there  is  a  one-parameter  family  of  unitary  operators 
U(t)   such  that 

(1.11)  U(t)"^L(t)U(t) 

is  independent  of  t.   This  fact  can  be  expressed  by 
setting  the   t   derivative  of   (1.11)  equal  to  zero: 

(1.12)  -U'^U^ULU  +  U"^L^U  +  U'-'-LU^  =  0 

A  one  parameter  family  of  unitary  operators  satisfies  a 
differential  equation  of  the  form 

(1.13)  U^  =  BU 

where  B(t)   is  an  antiself ad joint  operator.   Conversely, 
every  solution  of  (1.13 )  with  B*  =  -B,   is  a  one-parameter 
family  of  unitary  operators.   Substituting  (1.13)  into  (1.12) 

6 


we  get,  after  multiplication  by   U  ,on  the  left,   U    on 


the  right 


-  BL  +  L^  +  LB  =  0 


which  is  the  same  as 

(1.14)  L^  =  BL  -  LB  =  [B,L] . 

If  u   satisfies  the  equation  u^  =  K(u)  then  L^   can  be 
expressed  in  terms  of  u,   and  all  that  remains  to  verify 
is  that  equation  (l.l4)  has  an  antisymmetric  solution  B. 
The  drawback  of  this  method  is  that  it  requires 
one  to  guess  correctly  the  relation  (I.IO)  between  the  function 
u  and  the  operator   L.   Now  Gardner,  Kruskal  and  Miura  have 
made  the  remarkable  discovery,  [1],    that  the  eigenvalues 
of  the  SchrOdinger  operator 

(1.15)  L  =  D^  +  u/6 

are  invariant  if  u  varies  according  to  the  KdV  equation. 
We  shall  presently  verify  this  fact  with  the  aid  of  the 
linear  operator  equation  (l.l4);  more  generally  we  shall  use 
equation  (l.l4)  to  find  a  class  of  differential  equations 
under  which  the  operators   (1.15)    are  unitarily  equivalent 
for  all   t. 

With  the  choice  (1.15)  the  operator  L^   reduces 
to  multiplication  by  u^   so  that  according  to  (l.l4)  we  have 
to  find  an  antisymmetric  operator   B  whose  commutator  with 
L  is  multiplication.   An  obvious  choice  is 
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(1.16)  B   =  D 
^     o  o 

Indeed  an  easy  calculation  gives 

o^  -'    6   X 
This  shows  that  if  u  varies  according  to  the  equation 

(1.17)  ^t=^^x' 

the  operators  (I.I5)  are  unitarily  equivalent.   Alas,  this 
is  a  trivial  fact  since  changing  u   according  to  (I.17) 
amounts  to  replacing  the  potential  u  by  a  translate  of 
u,   which  obviously  results  in  an  equivalent  operator.   For 
a  less  trivial  result  we  try  a  third  order  antisymmetric 
operator 

(1.16)^  B^  =  D^  +  bD  +  Db, 

the  coefficient   b   to  be  chosen.   A  brief  calculation  yields 
the  following  value  for  the  commutator: 

[B-,,D]  =7tuD^+^u  D+iu 
1'  ■'    2  X     2  XX    6   XXX 

-  4b  D^  -  4b   D  -  b     +  4  bu 
X       XX     XXX   5   X 

Clearly  to  eliminate  all  but  the  zero  order  terms  we  have 
to  choose 

b  =  ^u 

With  this  choice   [B-,,L]   is  multiplication  by 
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-prr-[  U     +UU   ] 

24'-  XXX    x-' 

Multiplying   B-,   by   -24   we  verify  that 

[B^,D]  -  K(u) 

where   K(u)  =  u.   Is  the   KdV  equation  (1.1)1 

Clearly  this  process  can  be  generalized;  we  could 
choose  B   as  a  skew  symmetric  differential  operator 
of  any  odd  order   2q  +  1: 

B   =  D^q+l  +  ±   b  .D^J-l  +  D^J'-l  b  . 

q        1  J  J 

Since   [B  ,L]   is  symmetric  the  requirement  that  it  be  of 
degree  zero  imposes   q   conditions;  these  uniquely  determine 
the   q  coefficients   b .,   and  the  zero  order  term  of 
[B  ,L]  =  K   determines  a  higher  order  KdV  equation 

^t  =  Kq(^) 

which  shares  with  the   KdV  equation  the  property  that  the 
eigenvalues  of  the  Schrttdlnger  equation  with  u   as  potential 
are  its  integral. 

Gardner  has  discovered  an  interesting  relation  between 
the  higher  order   KdV  equations  and  the  explicit  sequence 
of  invariants  mentioned  earlier;  this  relation  will  be  describe 
in  section  2. 

The  process  described  above  can  be  generalized: 

Personal  communication 


Theorem  1.1;   Suppose  that   L   is  a  selfadjoint 
operator  depending  on   u   in  the  following  fashion; 

Lu  =  ^o  +  \ 


where   L    is  independent  of  u   and   M  depends  linearly 
on  u .    Suppose  that  there  exists  an  antisymmetric  operator 


B  =  B    such  that 
u  


^^'\^    =  \{u) 


Then  the  eigenvalues  of  L   are  integrals  of 

u,  =  K{u) . 

As  example  of  this  procedure  we  take  u  to  be  a 
p  X  p   syTumetric  matrix  variable  and  take   L  to  be  the  matrix 
operator   L  =  D   +  u/6.   If  we  choose   B  to  be  a  third 
order  matrix  operator  we  obtain  the  matrix  KdV  equation 

U,   +  -;=:(UU  +U  U)   +  U      =0. 

t    2^   XX      XXX 

Other  choices  for  the  operator   L   should  lead  to  other 
classes  of  equations. 

Having  shown  that  the  eigenvalues   A-,  (u)  ^  .  .  . ,  Ap.(u) 
SchrtJdinger  operator  (I.I5)  are  integrals  for  the  KdV 
equation  one  asks  how  these  integrals  are  related  to  the 
eigenspeeds   c-,  (u) ,  .  .  . ,  c^(u)   which  appear  in  the  asymptotic 
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description  (1.9).   Gardner  and  Kruskal  have  foiind  the  answer: 

(1.18)  c  .(u)  =  4a  .(u)  . 

J        J 

We  give  here  a  derivation  of  this  result  which  makes  use  of 
a  general  relation  involving  integrals  of  nonlinear  equations. 

This  relation  is  an  extension  of  the  following  well-known 
fact  about  quadratic  integrals   Q   for  linear  equations:  That 
the  bilinear  functional   Q(u,v)   is  independent  of   t   for 
any  pair  of  solutions   u  and   v.   This  result  can  be  deduced 
from  the  invariance  of  Q  for  the  superimposed  solutions 
u  +  V  and  u  -  v.   The  corresponding  result  for  integrals 
of  a  nonlinear  equation 

(1.19)  u^  =  K(u) 

is  derived  by  considering  one-parameter  families  u   of 
solutions;  these  can  be  constructed  e.g.  by  making  the 
initial  value  of  u  (t)   a  function  of  e: 

u^(0)  =  Uq  +  ef 

We  assume  that  the  nonlinear  operator  K  depends 
differentiably  on  u,   i.e.  that 


(1.20)  ^  K(u+ev) 


s=0  =  ^(-^^ 


exists  and  is  a  linear  function  of  v.   We  call  the  linear 
operator  V(u)   the  variation  of  K.   Differentiating  the 
equation 

11 


-rr-   u   =  K(u  ) 
dt   e      e 


we  obtain  the  variational  equation 
(1.21)  v^  =  V(u)v 

for  the  quantity 


du 
(1.22)   ,         v=^ 


e-0 


Let   I(u)   be  an  integral  of  the  equation  under  consid- 
eration;  we  assume  that   I(u)   is  dif ferentiable  in  the 
Frechet  sense,  i.e.  that 

-3—  I(U+£V)      ^ 

de   ^      e=0 

exists  and  is  a  linear  functional  of  v.   This  linear  functional 
can  be  represented  as 

(1.23)        ^  I(u+ev)  -  (G(u),v) 

where   (   ,   )   is  some  convenient  bilinear  functional; 
G(u)   is  called  the  gradient  of  I. 

Let  u  (t)   be  the  one-parameter  family  of  solutions 
considered  before;  then 

I(u^(t)) 

is  independent  of  t,      for  every  value  of   e.   Differentiating 
with  respect  to   e   we  obtain  that   (G(u),v)   is  independent 
of  t.   We  formulate  this  result  as 
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Lemma  1.1;   Let   u(t)   be  any  solution  of  equation 
(1.19),   v(t)   any  solution  of  the  corresponding  linear 
variational  equation  (1,21).   Let   I   be  any  integral  for 
(1.19),   G   its  gradient;  then 

(1.24)  (G(u(t)),  v(t)) 

is  independent  of  t. 

We  suppose  now  that  equation  (1.19)  is  translation 
invariant  and  that  it  has  a  solitary  wave  solution   s(x-ct). 
We  also  assume  that  the  bilinear  functional   (   ,   )   is 
symmetric  and  also  translation  invariant. 

Let  V  be  any  solution  of  (1.21);   according  to  (1.24) 

(G(s(x-ct),  v(x,t)) 

is  independent  of   t.   Using  the  translation  independence 
of   (   ,   )   and  introducing  the  abbreviation 

(1.25)  v(x+ct,t)  =  w(x,t) 

we  obtain  that 

(G(s(x)),  w(x,t)) 

is  equal  to  the  above  quantity  and  therefore  also  independent 
of  t.   Differentiating  with  respect  to  t  we  obtain  that 

(1.26)  (G(s),w^)  =  0. 
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From  (1.25) 

using  equation  (1.21)  and  the  fact  that  V  commutes  with 
translation  we  get 

w^  =  [cD+V(s)]w, 

where  D  denotes   S/Sx.   Substituting  this  into  (1.26)  we 
get 

(1.27)  (G(s),[cD+V(s)]w)  =  0. 

Denote  the  adjoint   of  V  by  V*    ;    since   (   ,   )   is 
translation  invariant,   D*  =  -  D,   so  we  get  from  (1.27)  that 

([-cD+V*(s)]G(s),w)  =  0. 

The  value  of  w  at  any  particular  time,  say  t  =  0  can  be 
prescribed  arbitrarily;  therefore  it  follows  from  the  above 
relation  that  in  fact 

(1.28)  [-cD+V*(s)]G(s)  =  0. 

Next  we  make  our  final  assumption  about  the  differential 
equation  (I.19):  it  is  energy  preserving,  i.e.  that  for 
solutions   u  of  (1.19) 

(u(t),u(t)) 

is  independent  of   t.   Differentiating  with  respect  to   t 
and  using  (I.19)  we  see  that'  this  amounts  to  requiring  that 
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for  all   u 

0  =  2(u,u^) =  2(u,K(u) ) 

The  initial  value  of  u   is  arbitrary.   Putting  u 

in  place  of  u  we  get  after  differentiating  with  respect 

to   e   that 

(v,K(u))  +  (u,V(u)v)  =  0. 

Using  the  adjoint  of  V  this  can  be  written  as 
(v,K(u)+V*(u)u)  =0. 

Since   v  is  arbitrary  this  implies  that 

(1.29)       K{u)  +  V*(u)u  =  0. 

We  turn  now  to  the  equation  for  solitary  waves: 

0  =  c  s   +  K(s) 
Expressing  K  from  (1.29)  we  can  rewrite  this  as 

[cD-V*(s)]s  =  0; 

i.e.  the  solitary  wave   s   belongs  to  the  nullspace  of  the 
linear  operator   cD-V*(s). 

We  assume  now  that  the  only  function  annihilated  by 
cD-V*(s)   and  vanishing  at  +  oo   is  a  multiple  of   s. 
From  (1.28)  we  see  that   G(s)   is  annihilated  by   cD-V*(s); 
therefore  if  G(s)   vanishes  at   +  oo   it  follows  that 
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Lemma  1.2;   Let   G  be  the  gradient  of  any  Integral 

of  (1.2);  then,  under  the  conditions  stated  before, 

V 

(1.30)  G(s)  =  Ks 

i.e.   that  every  solitary  wave  is  an  eigenfunction  of  the 
gradient  of  an  integral. 

It  is  easy  to  verify  that  the  KdV  equation  is 
energy  preserving  with  respect  to  the   Lp   scalar  product. 
In  that  case   K(u)  =  -uu   -u    ,   so  that 

V(u)v  =  -(uv)   -V 

^      '  ^     X    XXX 

and  thus 

V{u)  -  uD  +  eP. 

It  is  not  hard  to  show  that  every  function  annihilated  by 

(C-s)D-D^ 

and  vanishing  at   +  00   is  a  constant  multiple  of   s. 

We  apply  now  the  foregoing  to  the  integral   I(u)  =  ?\(u), 
where   A   is  an  eigenvalue  of  L  =  D  +  u/6: 

(1.31)  Lw  =  Tvw. 

To  compute  the  gradient  of   A  we  replace   u   by  u  +  ev 
and  differentiate  with  respect  to  e.   Denoting  -3—  by  a 
dot  and  using  the  fact  that   L  =  u/6  =  v/6  we  get 

Lw  +  vw/6  =  Aw  +  Aw. 

We  take  the   Lp   scalar  product  with  w;   using  the  fact  that 
L   is  symmetric  and  that  (I.31)  is  satisfied  we  get  rid  of 
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w   and  obtain 

(vw/6,w)  =  7\(w,w). 

Assuming  that  w  is  normalized  so  that   (w,w)  =   1 
we  get  that 

A  =  (vw/6,w)  -    /vw  /6  dx  =  (^  ,  v)  . 

*    d 
Since   A  =  —  I  =  (G(u),v),   this  shows  that  the  gradient 

de 

G   of   A   is  given  by 

2    ' 

(1.32)  G(u)  =  ^   . 

Since  eigenfunct ions  w  vanish  at   +  c©   the  hypothesis 
preceding  (I.30)  is  fulfilled;  therefore  we  conclude  from 
(1.30)  that   G(s)  is  a  constant  multiple  of  s: 

G(s)  =^=    J(s, 

i.e.  that  the  elgenfunction  w  of 

L  =  D^  +  s/6 
is 

(1.33)  w  =  const  s  ' 

This  relation  is  easily  verified  by  an  explicit  calculation; 
Taking  the  constant  to  be   1  we  get 

1     -1/2 
^x  =  2  ^x^ 


1  ^   ^-1/2    1  ^2^-3/2 

x^ 


^xx  =  2  ^xx^      -  W  ^-'^ 
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Using  relations  (1.6)  and  (1.7)  for  the  solitary  wave  we  get 

1/^    s^x  -1/2,  1/S^      2.       I 

Substituting  this  into  the  eigenvalue  equation 

T         ,  sw 
Lw  =  w   +  -7— 

XX     o 

we  get  after  a  brief  calculation  that 
r-,    ^h  \  T     T  1/2    c   1/2 

(1.3^)      Lw  =  Ls  ^   =  -q-  s  -^  . 

This  proves  that 

c(s)  =  4a(s) . 

Let  u  be  any  solution  of  the   KdV  equation  which 
contains  a  solitary  wave  travelling  with  speed  c,   i.e. 
such  that  given  any  positive  e   and  X,   there  exist  T 
such  that 

(1.35)  |u(x+cT,T)-s(x-e) I  <  e 

for  all   |x|  <  X.   We  claim  that  then  the  operator 

2 
Lpp  =  D  +  u(T)   has  c/k      for  an  approximate  eigenvalue  and 

w^(x)  =  s^/^(x-cT-9) 

as  approximate  eigenfunction,  in  the  sense  that 

(1.36)  I|LmW   -  7^  w||    <   6||w|| 
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where      5      tends   to    zero   as    e        — >     O'    and      X  — >■     cd  .      To 

see   this  we  use   the   fact   that   according   to    (I.35)    in  the 
interval 

(1.37)  cT-X<X<cT   +  X 

u  differs  by  e   from  s(x-cT-0);   therefore  using  (1.3^) 
we  conclude  that 

(1-38)        |L^w^  "  T  ^t'  ^  ^^t 

in  the  interval  (1.37).   On  the  other  hand  it  follows  from 
formula  (1.8)  that  outside  of  the  interval  (I.37)   w„  and 
its  second  derivative  are  bounded  by  exp( -const |x-x | ) . 
Denoting  by  M  the  supremum  u(x,t)   it  follows  that 

(1.38)'       |L^w^  -  ^   w^Ik  (const  +  M)e-^°''^^l^-^l 

outside  the  interval  (I.37).   Combining  (I.38)  and  (I.38)' 
we  deduce  (1 .36)  . 

According  to  spectral  theory  inequality  (I.36)  implies 
that   c/4   lies  within   6   of  a  point  of  the  spectrum  of  L„; 
since  we  have  seen  earlier  that  the  spectrum  of  L^  is 
independent  of  T,   it  follows  that   c/4   is  an  eigenvalue 
of  L. 

This  proves  one  half  of  (I.18);  the  second  half  - 
that  if  A   is  an  eigenvalue  of  L   then  4a   is  an  eigenspeed 
of  u  -  has  not  yet  been  demonstrated. 

In  section  2  we  shall  demonstrate  the  uniform  boundedness 
of  solutions  of  KdV  equation. 
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2.   In   [9]   Kruskal  and  Zabusky  have  studied  the 
Interaction  of  solitary  waves;  In  particular  they  posed  the 
following  problem:  \ 

■'   Let   d  be  a  solution  of  KdV  which  for   t   large 
negative  represents  two  solitary  waves  travelling  with 
speed   C-,  ,   respectively   Cp,  approaching  each  other; 
what  Is  the  asymptotic  behavior  of   d(x,t)   for  large 
positive  values  of  t? 

They  solved  this  problem  by  computing   d;   for   -T 
large  negative  set 

(2.1)         d(x,-T)  =  s(x,c-j^)  +  sCx-XjCg) 

where  the  separation  distance  X  was   chosen  so  large 
that  the  two  solitary  waves  overlap  only  by  a  negligible 
amount;  since  solitary  waves  die  down  exponentially,  even 
a  moderate  value  of  X  accomplishes  this.   The  speeds  were 
chosen  so  that   c-,  >  c^,   and  X  was  taken  to  be  positive 
so  that  at  t  =  -  T  the  faster  wave  lies  to  the  left  of 
the  slower  one . 

If  the  equation  governing  the  motion  were  linear,  the 
solitary  waves  would  not  interact  at  all  and  so  after  the 
elapse  of   2X/(c-,-Cp)   time  the  relative  position  of  the 
two  would  be  merely  interchanged.   Numerical  calculation 
of  the  solution  of  the  KdV  equation  with  initial  values 
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(2.1)  showed  that  for   S  >  2X/c-^^-C2 

d(x,-T+S)  =  s(x-c-^S-0-j_,c-j_)  +  s(x-X-C2S-02,C2), 

except  for  deviations  that  could  be  accounted  for  by- 
truncation  error,   i.e.  the  same  as  would  be  in  a  linear 
case  except  for  phase  shifts   9-,   and   9p. 

The  actual  process  of  interaction,  i.e.,  the  behavior 
of   d(x,t)    around  the  time   -  T  +  X/c-, -Cp   is  far  from 
being  a  mere  superposition.   In  fact  Kruskal  and  Zabusky 
observed  that  in  cases  when   c,  »   Cp,  i.e.  when  the  first 
wave  was  much  higher  (and  therefore  faster)  than  the  second 
one,  the  big  wave  swallows   up  the  small  one  during  the 
interaction,  and  reemits  it  later.   In  cases  where   c,   and 
Cp  were  comparable  the  two  waves  interact  as  follows:   As 
soon  as  the  big  wave  comes  reasonably  close  to  the  smaller 
one  in  front  of  it  the  big  wave  begins  to  shrink  and  the 
smaller  one  begins  to  grow,  until  the  two  waves  interchange 
their  roles;  thereafter  they  separate. 

In  this  section  we  present  a  rigorous  proof  that  the 
KdV  equation  indeed  has  solutions  which  behave  like  two 
solitary  waves  approaching,  interacting,  then  separating; 
furthermore  we  give  a  precise  estimate  for  the  ratio  of  speeds 
which  lead  to  the  two  different  kinds  of  interaction  described 
above.  It  turns  out  that  there  is  yet  a  third  manner,  inter- 
mediate between  the  other  two,  in  which  the  interaction  can 


21 


take  place. 

We  shall  call  these  special  solutions  double  waves, 
and  denote  them  by  d(x,  t;  c-,  ,  Cp)  .   Note  that  for  fixed  c, 
and   Cp  there  is  a  two-parameter  family  of  double  waves, 
the  parameters  being  the  phase  of  each  solitary  wave  at 
t  =  -  CO  . 

In   [8]    Kruskal  and  Zabusky  derived  an  ordinary 
differential  equation  with  respect  to   x  which   d(x,t;c-.Cp) 
satisfies  for  each  fixed   t.   The  argument  in   [8]    is 
formal;   here  we  present  an  (almost)  rigorous  derivation. 

Our  starting  point  is  lemma  1.1,  the  constancy  of 
the  functional  (1.24)  for  pairs  of  solutions   u,  v   of 
the  nonlinear  equation  and  of  the  variational  equation 
respectively.   Taking  u  to  be  d  we  get  that  for  any 
solution  V  of  (1.21), 

(2.2)  (G(d),v) 

where   G   is  the  gradient  of  some  integral,  is  independent 
of  t.   By  definition,  for   t   large  negative 

(2.3)  d(x,t;c-|^C2)  =  s  (x-c^t-Q-j^,  c-j^) 

+  s(x-C2t-92,C2) 
+  error (t ) . 

The  error  term  in  (2.5)  is  due  to  the  interaction  of  the 
tails  of  the  solitary  waves  into  which  the  double  wave   d 
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decomposes.   Since  the  tails  of  solitary  waves  decay 
exponentially,  and  since  the  separation  of  the  two  solitary 
waves  is  proportional  to   t,  we  expect  the  error  term  in 

(2.3)  to  decay  exponentially  in   t.   Suppose  that  the 
gradient   G   is  a  .local  operator;  it  follows  from  (2.3) 
that 

(2.4)  G(d)  =  G(s^)  +  GCs^)  +  error, 

where  the  error  in  (2.4)  also  tends  to  zero  exponentially 
as   t  — ^  -  00  . 

Suppose  that   G   satisfies 

(2.5)  G(s^)  -  0,  G(S2)  =  0; 

then  it  follows  from  (2.4)  that  V(d)   tends  to  zero 
exponentially  as   t  — >  -  oo  . 

¥e  turn  now  to  the  fiinctions   v;  since  these  satisfy 
a  linear  equation  it  can  be  shown,  by  an  argument  similar 
to  the  one  which  led  to  inequality  (1.4),  that  any  v 
increases  at  most  exponentially.   In  fact  since  exponen- 
tially increasing  solutions  of  the  variational  equation 
usually  indicate  an  instability,  and  since  on  the  other 
hand  numerical  evidence  Indicates  that  double  waves  are 
stable,  it  is  reasonable  to  expect  that  all  solutions   v 
of  the  variational  equation  grows  at  a  rate  slower  than 
exponential . 

We  return  now  to  the  functional  (2.2);  we  have 
shown  that  as   t  — >  -  go   the  first  factor   G(d)   tends 
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to  zero  exponentially,  and  that  the  second  factor   v   tends 
to  infinity  slower  than  exponentially.   It  follows  then 
that   (G(d),v)   tends  to  zero  at   t  — >  -  oo  ;   on  the  other 
hand  according  to  (2.2)  this  functional  is  independent  of 
t.    Therefore  it  follows  that   (G(d),v)   is  zero. for  all   t, 
At  any  particular  time  t   the  initial  values  of  v  may 
be  prescribed  arbitrarily;  therefore  it  follows  that 

(2.6)  G(d)  =  0 

at  time   t  ,   that  is  for  anytime. 

There  remains  the  task  of  constructing  an  integral  I 
whose  gradient   G  is  local  and  satisfies  (2.5),  i.e. 
annihilates  both  solitary  waves   s.   and   Sp.   Here  we 
rely  on  lemma  1.2,  relation  (I.30),  according  to  which 
solitary  waves  are  eigenfunctions  of  the  gradient  of 
every  integral: 

G(s^)  =X(G)s-L   ,   G(S2)  =><{G)S2 

It  follows  that,  given  three  independent  integrals, 
an  appropriate  linear  combination  of  them  will  have  a 
gradient  which  annihilates  any  two  given  solitary  waves. 

We  turn  now  to  the  task  of   finding  three  independent 
integrals  whose  gradients  are  local  operators .   We  have 
already  noted  in  section  1  that  the  energy 

(2.7)-L        I-l(u)  =  Ti  u^dx  -  |(u,u) 
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is  an  Integral  for  the   KdV  equation.   Another  integral 
was  found  by  Whitham,  [ 7] : 

{2.7)2  ^2^^^  ^I^J"^^    ~   ^x^^^""- 

A  third  integral  was  discovered  by  Kruskal  and  Zabusky: 

(2.7)3        ^jt"'  =/'?"'  -  J^x'  -  f  "xx'>^- 

Subsequently  Kruskal  and  Zabusky  found  two  more  explicit 
integrals  and  Miura  four  more;  an  infinite  sequence"*"  of 
them  was  constructed  in  [5]. 

For  the  discussion  of  the  double  wave  we  need  only 
the  first  three  integrals.   The  gradients  of  these  are 

(2.8)^        G^(u)  =  u, 

(2.8)2        ^2^^^  =  u^  +  2u^^, 

(2.8)3        ^3(-)  --'  ^  V-6-xx-f  -xxxx- 

Note  that  these  are  indeed  local  operators.   The  proof 
that   I^,  1^,    I   are  integrals  of  KdV  consists  in 

1  th 

Denote  by  I^  the  n   integral  of  this  sequence  and 

by  G   its  gradient.   The  relation  discovered  by 

th 
Gardner  between   I^   and  the  n   generalized  KdV 

operator   K^   described  in  section  1  is  as  follows: 

^n  =  D^n- 
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verifying  that  the  product  of  G  (u)   with  K(u)  =  uu  +  u 
is  a  perfect   x   derivative.   Indeed  an  explicit  calculation 
gives  \ 

{2.9)-L     G-l(u)K(u)  =3  (uu^^  _  1  u^^  +  i  u5)^  =  H-L^^^x 

(2.9)2     G2{u)K(u)  =  (u^^  +  u^u^  +  ^   u^)^  =  H^Cu)^ 

(2.9),     G^(u)K{u)  =  (|-u^  +^u   UU  + 

+  ^U    u-l-u    U    +4u   U  +  U   U'^+^) 
5  XX     5  XX  X    2  X      XX     5   x 


=  H^(u) 
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We  know  that  solitary  waves  are  eigenfunctions  of 
G  ;   a  calculation  gives  the  eigenvalues  }1{G    )    as  follows: 

(2.10)-,  G  (s)  =  s 

1  1^ 

{2.10)2  ^2^^^  "  ^^^ 

(2.10)^  G^(s)  =  ^  c^s 

3  3      5 

We  form  now  the  linear  combination 

(2.11)  I  =  I^  +  AI2  +  BI-^ 

whose  gradient  is 


(2.12)  G  =  G,  +  AG2  +  BG^ 
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In  view  of  (2.10) 

G(s)  =  (^  c^  +  2Ac  +  B)s. 

Thus   G  annihilates   s-^   and   s^   if   c^   and  c        satisfy 
the  equation 

(2.13)      S(c)  =  ^  c^  +  2Ac  +  B  =  0. 

In  view  of  the  relation  between  coefficients  and  roots  of 
a  quadratic  equation  this  means 

(2.1^)^  A  =:  -  ^(c^+c^). 

In  view  of  (2.6)  we  conclude:  Let   d  be  a  double 
wave  with  speeds   c-j^  and  Cp,   and  define  the  constants 
A,  B  by  (2.14);  then  for  each  t,   d   satisfies 

(2.15)      G(d)  =  G^(d)  +  AG2(d)  +  m^{d.)    =   0 

(2.15)  is  a  nonlinear  ordinary  differential  equation 
of  fourth  order  ;  therefore  its  solutions  form  a  4 -parameter 
family  of  fimctions;  on  the  other  hand  the  double  waves  form 
a  2-parameter  family.   We  shall  accordingly  deduce  from 
(2.15)  a  second  order  equation  satisfied  by  all  double  waves. 
To  do  this  we  make  use  of  the  earlier  observation,  see  (2.9), 

This  equation  appears  in  [8]. 
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that  the  invariance  of  I(u)   for  solutions  of  KdV  is 
equivalent  with  the  fact  that   G(u)K(u)   is  a  perfect  x 
derivative.   Multiplying  (2.15)  by  K(u)   and  using  (2.9) 
we  get  that  ^ 


(H^  +  AH2  +  BH-^)^  =  0, 


Integrating  this  and  using  the  fact  that   d  and  all  its 
derivatives  are  zero  at   x  =  +  00   we  deduce 

(2.16)        H^(d)  +  AH2(d)  +  BH^(d)  =  0. 

Next  we  make  use  of  the  translation  invariance  of 
the  integrals   I  under  consideration,  i.e.   that   I(u  ) 
is  independent  of  e,      where  u   is  the  translate  of  u 
in  the  x  direction  by  e .   Differentiating  with  respect 
to  e   we  get 

0=^I(u^)|^^0-/^(-)^x  ^^- 

This  implies  that   G(u)u   is  a  perfect   x   derivative. 
Indeed  an  explicit  calculation  gives 

(2.17)-L    ^l^^^'^x  =  (|  ^^^x  =  '^l^^'^x 
(2.17)2    G2(u)u^=  (u/+|^)^  =  J^{u)^ 

h 

(2.17)^     G,  (u)u  =    (^  u  u     -  I-  u   +  3uu   +  7^) 
^    '  '3     3    X   ^5   X  XXX   5   XX   -^  X    4  'x 

Thus  multiplying  (2.15 )  by  u   we  get,  after  inte- 
gration  and  using  the  fact  that   d  and  its  derivatives 
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vanish  at   oo  ,   that 

(2.18)  J^(d)  +  AJgCd)  +  BJ^(d)  =  0 

Both  (2.16)  and  (2.l8)  are  third  order  differential 

equations.   Expressing   d^^^   from  (2.l8),  substituting 

2 
into  (2.16)  and  multiplying  by   d^   we  get  an  equation 

of  second  order  and  fourth  degree  in   d    which  we  write 

symbolically  as 


(2.19)  Qt^xx'^x'*^^  =  °- 


Prom  this  equation   d    can  be  expressed  as  a  4-valued 
function  of   d   and   d  >   since   x   does  not  appear 
explicitly  in  these  equations  this  second  order  equation 
is  equivalent  to  a  first  order  autonomous  system  of 
equations  for   d  and   d  .   By  studying  carefully  the 
geometry  of  all  four  branches  of  the  corresponding  vector - 
field  one  can  show  that  (2.19)  indeed  has  solutions  whj.ch 
tend  to  zero  as   x  -^  +  00   and  has  the  shape  of  a  double 
wave,  i.e.  has  two  maxima  and  one  minumum.   We  shall  not 
present  the  details  because 

a)  We  didn't  carry  them  out  completely 

b)  The  formulas  are  horribly  complicated 

c)  An  explicit  formula  for  double  waves  (indeed 
N-tuple  waves)  was  derived  recently  in   [l]. 

We  show  now  how  to  use  equation  (2.19)  to  study  the 
time  history  of  double  waves .   For  this  purpose  we  study 
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the  maximum  value  of   d(x,t)  with  respect  to  x,   or 
rather  the  relative  maxima  of  d  as  functions  of  time. 
Denote  by  m  =  m(t)   the  value  of  a  relative  maximum  of 
d(x,t);   at  the  point   y  where  the  relative  maximum  occurs 


(2.20)  ^x  "  ° 


It  follows  from  (2.20)  by  the  implicit  function  theorem 
that  if   d   <  0  at   y,   y  is  a  differentiable  function 
of  t;   then  m  =  d(y,t)   also  is  differentiable  and 
satisfies 

"^t  =  ^^x^t  +  ^t  ^  ^t    • 

Since   d   satisfies  the   KdV  equation,  we  have,  in  view 
of  (2.20) 


(2.21)        m^  =  -  K(d)  =  -  d^^^. 

We  proceed  now  to  determine   d     at  a  local  maximum. 
At  a  point  where   d  =  0  equation  (2.l6)  simplifies  consid- 
erably;  using  formulas  (2.9)  and  denoting  the  value  of  d 
by  m  we  get 

9,2   ,6,2^^^  J>   ^rrp 


(2.22)        F-d    +^dm  +  d  m-^    +  -^ 
^     '        5  xxx   5   XX     XX     5 


+  A(d  ^  +  m^d    +  i  m  )  +  B(md^^  +  ^  m^)=0, 
^  XX        XX    4     '      ^   XX    3 
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Similarly,  at  a  point  where   d  =0   equation  (2.l8) 
becomes 

(2.23)  d^2  ^  l^m  ^  ^  m  ^  ^  m^^  ^  p^^^ _ 

From  (2.23)  we  deduce  that 

1/2 

(2.24)  d^^  -  -  P(m); 

the  negative  sign  is  chosen  on  account  of  the  second 
derivative  being  nonposltive  at  a  local  maximum  point 
Substituting  (2.24)  into  (2.22)  we  get 


R(m)    =   -   a(m)    +  b(m)P^/^(m), 


(2.25) 

where 

(2.26) 

with 

(2.27) 

and 

(2.27) 

f    \         2     ^/    X     ,    m^,    5AP(m)     ,    5Am     ^  5B     3 
(m)    =  r=r  mP(m)    +  q   +  ^q    ^        +  ^-rz—  +  ^  m^ 


^  <A{mj    =  J  iiir^iu    +   9-  -1-  -9     1-     3fc      -^   27 


b(m)    =  ^  Tcr    +  ^  m     +  4  Bm 


Combining  (2.21)  and  (2.25)  we  have 

(2.28)  ^t  =  ±  R("i)^'^^; 

the  sign  to  be  taken  as  positive  when   m  increases, 
negative  when  m  decreases.   Thus   m  as  function   t 
is  governed  by  equation  (2.28). 

To  study  the  behavior  of  m  we  have  to  know  something 
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about  the  function  R(m).   We  turn  now  to  this  task;  we 
start  with  the  observation  that  the  solitary  waves   s, 
and   Sp  themselves  satisfy  the  double  wave  equation, 
and  therefore  their  maxima  satisfy  equation  (2.28). 
On  the  other  hand  the  maximum  of   s, ,   respectively 
Sp,   is  independent  of  t,   being  equal  to 

(2.29)  m-,  =  3c-,,  respectively  mp  =  3Cp. 

This  shows  that  the  constants   m-,  ,   mp   satisfy  (2.28) 
which  implies  that  R  vanishes  at  m,   and  mp.  Actually 
more  than  this  is  true: 

Lemma  2 . 1   R(m)   has  a  double  zero  at  m,   and  at 

d^R 
mp,   and  — t^     is  positive  at  these  points . 

dm 

Proof;   The  polynomials  P  and  S,   defined  by 

(2.23),  respectively  (2.13),  satisfy 

(2.30)  P(m)  =^m2s(m/3)  +  j^   . 


Since  S(c)  =  0  for   c  =  c,   or  Cp,   it  follows  by 

(2.29)  and  (2.30)  that  P(m)  =  m  /36  for  m  =  m^  or 

mp,   and  so  for  these  values 

1/2     2 
P(m)  =  m  /6. 

Using  this  relation  and  expressions  (2.l4)  for  A  and 
B  we  can  express   dR/dm  and   dT^/dm   at   m,   and  mp 
as  polynomials  in  m,   and  mp.   It  is  then  a  matter  of 
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explicit  calculation,  omitted  here,  to  show  that  the 
former  are  zero  and  the  latter  positive.   This  completes 
the  proof  of  lemma  2.1.   It  follows  that  R(m)   is 
positive  in  some  neighborhood  of  m,   and  m_;   we  turn 
now  to  the  behavior  of  R(m)   in  the  whole  interval 
(m-,,mp).   First  of  all  we  investigate  the  sign  of  P{m); 
since   P(m)   is  ^  m   times  the  quadratic  polynomial 

(2.31)  ^  +1  m  +  I  , 

it  will  be  positive  if  the  discriminant  of  (2.31)  is: 

2 
B   A^ 
discr  =  -p   -  -q 

Using  formulas  (2.14)  for  A  and   B  we  get 

9  2      2 

discr  =  -^(^c-j^Cg-  ^i   "  ^2  ' 

This  quadratic  form  is  positive  if  and  only  if 

(2.32)  -^  1^^^  =  2.62 

Thus  we  have  proved 

Lemma  2.2:   If   c.^   and   c^   satisfy  (2.32),   P(m) 

is  positive  for  all  real  values  of  m.   If  (2.32)  is 

violated,   P(m)   is  negative  in  the  interval   (n-|^,n2), 

(2.33)  "12"  f["^i+"^2i  (m^^+in2^"^'^'"^2^    ^ 

It  is  easy  to  verify  that  the  interval   (n  ,n2) 
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lies  inside   (iii^,mp). 

We  turn  now  to  R(m);   adding   a(m)   to  both  sides 
of  (2.26)  and  squaring  we  get  ] 

(2.34)  R^  +  2Ra  =  b^P  -  a^  =  T(m)m  . 

Clearly   T   is  a  polynomial  in   m  of  degree   6.   It  ■ 
follows  from  (2,34)  that   T   vanishes  wherever  R   does. 
Since  R   has  a  double  zero  at  each  of  the  points   m-, 
and  mp,   so  does   T;   a  large  amount  of  highschool 
algebra  yields  the  remaining  quadratic  factor  of  T  and, 
just  as  in  highschool,  that  quadratic  polynomial  can  be 
factored.   The  resulting  factored  form  of  T   is 

2       2 

(2.35)  T(m)  =  const(m-m-j^)  (m-m^)  (m-m-j^-m^)  (m+m-j^-m^) 


Lemma  2.3;   a)   In  the  range 

c  /c   <  ^  +v^ 
12     2 


the  function  R(m)   is  positive  in   (m-,,m2) 

b)  In  the  range 

^  +\/5  <  c  /c   <  3 
2      r  2   ^ 

R   is  positive  in   (m^,n-j^)   and  in   (m^jn^), 
n, ,  np  are  defined  by  (2.33). 

c )  In  the  range 

3  <  c^/c^ 
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where 


R 


Is   positive    on      [m^,n^)      and   on      (m-j^jm-j^-m^) ,      and 


R{m-,  -iTip)    =    0, 


a) 


R 


b) 


R 


c) 


R 


Fig.    1 
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Proof;  As  remarked  earlier  it  follows  from  (2.3^) 

that  every  root  of  R   is  also  a  root  of  T.   It  follows 

from  (2.35)  the  only  possible  root  of  T   on   (m-,,mp) 

is   m  =  m-,  -  mp.   An  explicit  substitution  into  (2.23) 

and  (2.27),  using  relations  (2.l4),  gives 

1  2 

Plm^-mp)  =  ^  (m^-3m2), 

1  2 

(2.36)    a(m)  =  -^  (m-|^-2m2)  (m-j^-3m2), 

b(m)  =  -g-   (m-|^-2m2)  . 

From  the  first  of  these  relations  we  have  that 

-,  /p  -,  -7-  (3ii^p-iTi-|  )  in  cases  a)  and  b) 

P  /  (m^-mp)  =  -^  |m-|^-3m2|  = 

■7-  (m, -3mp)  in  case  c) 

Substituting  this  into  the  definition  (2.26)  of  R  we  get, 
after  a  short  calculation,  that 

^  0     in  cases  a)  and  b) 


R(m^-m2) 


0  in  case  c). 


as  asserted.   It  is  easy  to  verify  that  in  case  c), 

m-,  -  mp  is  greater  than  n  ;   this  completes  the  proof 

of  lemma  2.3 . 

Eventually  we  will  be  interested  not  only  in 
the  value  of  the  maximum  as  function  of  t   but  also  in 
its  location  y(t): 

d(y(t),t)  ^   m(t) . 

The  criterion  for  a  maximum  is 
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d^(y,t)  =  0; 

differentiating  this  with  respect  to   t   we  get 

XX  "^  t     xt 
or 

(2.27)  yt  =   -  T^    • 

XX 

Differentiating  the   KdV   equation  we  get 

(2.28)  d^^  =  -  d/  -  dd^^  -  d^^^^  . 

Using  equation  (2.15)  we  can  express   d      as  function 

of  d,  d  ,  d   .   and   d   ;  furthermore,  at  a  local  maximum 
'   x'   XX'         xxx' 

point  we  have   d  =  0,   and  d   ,  d     can  be  expressed 

X  XX     XXX 

as  in  (2.24)  and  (2.25)  as  functions  of  m.   Substituting 
these  expressions  into  (2.28)  and  (2.27)  we  get  an  expression 
of  the  form 

(2.29)  Yt  =  Y(m). 

If  we  know  m  as  function  of  t,    y(t)   can  be  determined 
by  integration. 

We  shall  not  calculate   Y(m)   explicitly;  it 
suffices  to  note  that  the  discussion  applies  to  the  solitary 
waves   S-,   and   Sp;   for  these  the  value  of  the  maximum  is 
m.,  ,   respectively  m^,      and  the  location  of  the  maximum 
moves  with  speed  c-j^  ""  J  ^i      ^^'^     ^2  "^  T  "^2   ^-espectively . 
Therefore  (2.29)  implies  that 
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(2.30)         Y(m^)  =  c-^  ,   Y(m2)  =  c^. 

We  are  now  in  a  position  to  prove  \ 

Theorem  31:   For  any  pair  of  speeds  c-,   and   Cp 

there  exists  a  double  wave,  i.e.  a  solution  d(x,t)   of 

the   KdV  equation  such  that 

(2.31)  d(x,t)  -  s(x-c^t-0i;c-j^)  -  s(x-C2t-0i;C2) 

tends  to  zero  uniformly  as   t  — >  +  oo  . 

Proof;   We  take  as  initial  value  of  d  a  solution 
(2.19);  according  to  the  existence  theorem  in  [6]  the 
KdV  equation  has  a  solution   d(x,t)   for  all  time  with 
these  initial  data,  and   d  and  all  its  derivatives  are 
zero  at  x  =  +  00  .   We  claim  that  for  all  values  of  t, 
xi(x,t)   satisfies  (2.19).   To  see  this  we  note  first  that 
since   d(x,0)   satisfies  (2.19),  it  also  satisfies  (2.15): 

(2.32)  G(d(0))  =  0. 

According  to  lemma  1.1,  equation  (1.24),  for  any  solution 
V  of  the  variational  equation 

(2.33)  (G(d),v) 

is  independent  of   t.   Relation  (2.32)  shows  that  at 
t  =  0,   (2.33)  is  zero;  therefore  (2.33)  is  zero  for  all 
time.   Since  the  value  of  v   can  be  prescribed  arbitrarily 
at  any  particular  time,  it  follows  that   G(d)   is  zero  at 
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each  time.   From  this,  and  the  fact  that   d   and  its 
derivatives  are  zero  at   x  =  Hh  oo   we  deduce  as  before 
that   d   satisfies  equation  (2.19)  at  each   t. 

We  turn  now  to  case  a)  and  claim: 

For  any  time   t,  d(x,t)  has  exactly  two  local  maxima. 

Proof;  In  case  a)  the  number  of  local  maxima  is 
independent  of  t  since  at  the  time  of  the  creation  or 
disappearance  of  a  local  maximum  d   =  0;   on  the  other 

XX 

hand,  according  to  (2.23)   d   =  P(m)   at  a  local  maximum, 
and  in  case  a)   P(m)   is  positive  for  all   m. 

The  time  history  of  each  local  maximum  is  governed 
by  equation  (2.28).   In  case  a)   R   is  positive  in   (m-,,mp) 
and  has  a  double  zero  at   m   and  at   m  ;   it  follows  from 
this  that  each  solution  of  (2.28)  whose  values   lie  in 
(m-,,mp)   goes  from  m-,   to  mp  as   t   goes  from  +  oo  to 
-  oo,   or  the  other  way,  depending  on  the  sign  in  (2.28). 
Furthermore  the  approach  of  m(t)   to   m-,   or   mp  as 
t  — ^+   CD   is  exponential.   It  follows  from  this  that  the 
function  Y(m)    occurring  in  (2.29),  (2.^0)  tends  to   c-,, 
respectively   Cp  exponentially  as   t  — >  _+  oo  .   Integrating 
(2.29)  we  conclude  therefore  that 


(2.3^)  y(t)  -  c^^2^ 


It  is  not  hard  to  show  that  the  values  of  any  solution 
of  (2.19)  which  is  zero  at   x  =  +  ^   ^^^  ^^    [m-^,m^)  . 
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tends  to  a  limiting  value  as   t  — >  +  °° • 

Next  we  show  that  as   m(t)  — ^  m-,   or  m^,      the 
shape  of  the  curve   d(x,t)   around  the  maximum  point  y(t) 
tends  to  the  shape  of  a  solitary  wave,  i.e.  that 


(2.35)       lim  d(x-y(t),t)  =  s^  ^{x) 
t— >  +00  ■* 


t->_ 

uniformly  on  bounded  x-intervals.   To  prove  this  we  merely 
note  that  for  each  fixed   t,  s-|^(x),  s^{x)      and   d(x,t) 
all  satisfy  the  fourth  order  equation  (2.I5).   Furthermore 
it  follows  from  relations  (2.20),  (2.24),  (2.30)  and  (2.25) 
that  the  Cauchy  data  of  d  at   x  =  y   for  a  fourth  order 
ordinary  differential  equation,  i.e.  the  values  of   d  and 
of  its  first  three  x-derivatives,  tend  as   t  — >•  +00   to 
the  Cauchy  data  of   s-,(x),   respectively   S2(x)   at   x  =  0. 
Relation  (2.35)  follows  then  from  the  continuous  dependence 
of  solutions  on  their  Cauchy  data. 

From  (2.35)  we  can  deduce  our  assertion  about  the 
number  of  maxima;  for  assume  that  these  were  more  than  two. 
Then  at  least  two  of  them  would  tend  to  the  same  limit, 
say  m-,,   as   t  — >  00  .   It  follows  from  (2.34)  that  the 
separation  of  these  two  locations  of  maxima  tends  to  a 
constant;  but  this  is  incompatible  with  relation  (2.35) 
which  says  that   d  looks  like   s-,   centered  around  either 
location . 

Likewise  it  is  impossible  for   d   to  have  only  one 
maximum;  for  then   d(x-y(t),t)   would  be  monotonically 
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decreasing  in  x  on  either  side  of  ■  y(t)  and  so  would 
tend  uniformly  to  one  solitary  wave  as  t  — >  -  od  ,  the 
other  as   t  —^   +00.   We  claim  that  for  all   t 

(2.36)  J  d^(x,t)  dx  <  e(X), 

|x-y(t) |>X 

where   e(X)   tends  to  zero  as   X  — ^  00  .   If  this  were  so 

we  could  deduce  that 


J^l    (x)dx 

(2.57)     11m     rd^(x-y(t),t)dx  = 

t— >  +  00  '-'  p     2 

/sg  (x)dx 


But  this  Is  a  contradiction.   For,  on  the  one  hand,  for 

/s.  dx  7^  /s^  dx  ; 

on  the  other  hand   /d  dx   Is  Invariant  under  both  x  and 
t   translation,  and  so  the  integral  on  the  left  side  of 
(2.37)  is  Independent  of  t;  but  then  the  limits  as 
t  — >  +  00   and   -  GO   cannot  have  different  values. 

To  prove  (2.36)   we  merely  note  that 
/d(x,t)dx  =  IqC^)   is  also  an  Integral  for  the   KdV 
equation,  from  which  (2.36)  follows  with 

e(X)  =  s(X), 

because  of  the  monotonlclty  of  d   on  either  side  of  y(t) 

Having  shown  that   d(x,t)   has  exactly  two  maxima 
it  follows  easily  by  the  arguments  already  presented  that 
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as   t  —^   +00,  d(x,t)   tends  uniformly  to  the  superposition 
of  two  solitary  waves,  each  travelling  at  its  own  speed. 
This  completes  the  proof  of  theorem  2.1  in  case  a). 

The  analysis  presented  above  shows  that  the  time 
history  of  the  two  maxima  is  as  follows:  as   t   goes  from 
-  oo   to  +  oo  the  height  of  the  larger  solitary  wave 
decreases  from  m-,   to  mp,   while  the  height  of  the  smaller 
one  increases  from  m^  to  m-,  .   Thus  in  this  case  the  two 
solitary  waves  interchange  their  roles  without  passing 
through  each  other,  as  observed  by  Kruskal  and  Zabusky 
in  their  calculations. 

We  turn  now  to  case  b);  here  R(m)   is  not  defined 
inside  (n, ,np),   which  means,  see  equation  (2.23),   that 
no  relative  maximum  of  d(x,t)   can  lie  in  that  interval. 
We  assert  that  the  absolute  maximum  lies  above  that 
interval  at  all  times.   For  suppose  on  the  contrary  that 
at  some  time  the  absolute  maximum  is  less  than  n^;   then, 
since  the  absolute  maximum  is  a  continuous  function  of  t, 
and  since  on  the  other  hand  the  value  of  the  absolute 
maximiun  cannot  cross   (n-,,np),   it  follows  that   d(x,t) 
does  not  exceed  np  for  any  value  of  x  and  t.   Let 
m(t)   denote  a  local  maximum  at  time   t.   As  we  saw 

earlier  m(t)   satisfies  the  differential  equation  (2.28): 

1/2 
m.  =  +  R    (m)  •   I"^  case  b)   those  solutions  of  this 

differential  equation  whose  values  lie  in   (mp,np)   behave 
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as  follows:  depending  on  the  sign  in  (2.28),   m(t)   tends 
with  increasing  (decreasing)   t   to  np: 

lim  m(t)  =  n„  . 

t->  t 
^  o 

Denote  as  before  the  location  of  the  maximum  by  y(t), 

and  set 

lim  y(t)  =  y  . 

t->  t         ° 
^     o 

It  follows  from  relations  (2.20)  and  (2.24)  that 


d  (y  ,t  )  =  0,   d   (y  ,t  )  =  0 
x^^^o'  o'  '    xx^'^o'  o 

while  it  follows  from  (2.25)  and  the  fact  that 
R(n2)  >   0  that 

d    (y  ,t  )  >  0. 
XXX  ^  -^  o '  o ' 

This  implies  that  the  function  d(x,t  )   does  not  have  a 
local  -  much  less  global  -  maximum  at   x  =  y  .   Since 
the  value  of   d   at   x  =  y   is   n^,   it  follows  that  the 
global  maximum  of   d(x,t  )   exceeds   n^;   this  contradicts 
our  previous  assumption,  and  so  proves  the  assertion. 

Denote  by  M(t)   the  absolute  maximum  of   d(x,t); 
denote  by  m(t)   that  local  maximum  which  at,  say,  t  =  0 
equals   M(0)   and  which  satisfies  the  differential  equation 
(2.28)  for  local  maxima.   The  function   m(t)   is  either 
increasing  or  decreasing;  suppose  it  is  increasing,  then 
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since  all  local  maxima  satisfy  (2.28),  it  follows  that 

m(t)   remains  the  absolute  maximum  for   t  >  0;   it  further 

follows  from  (2,28)   that   m(t)   tends  to  m,   as   t 

tends  to   +  oo  ,   and  that  for  some  finite  time   t 

lim  m( t )  =  n, . 

When  m  reaches  the  value   n-,  ,  it  ceases  to  be  a  local 


maximum,  so  that 


MCt^)  >  m(tQ). 


Denot 


e  by  t-i   the  infimum  of  those  values  of  t-,  for 


which  m(t)   is  the  absolute  maximum;  clearly   t   <  t  <  0. 
At  time   t-,   there  must  be  at  least  two  points  where  the 
absolute  maximum  is  assumed;  denote  by  n(t)   that  local 
maximum  which  is  equal  to  m(t)   at   t  =  t,   but  is  a 
decreasing  function  of  t;   n(t)   satisfies  equation  (2.28)_, 
and  it  follows  that  as   t  — >  -  co ,  n(t)   tends  to  m-,  . 

Denote  by  y(t)   the  location  at  time   t   of  the 
absolute  maximum.   As  before  we  deduce  that 

lim    d(x+y(t),t)  =  s-,(x), 
|t|-^  oo 

uniformly  on  bounded  x   interval. 

The  analysis  presented  in  the  previous  case,  when 
applied  to  this  situation,  shows  that  for   |t|   large  there 
can  be  at  most  one  additional  local  maximum;  suppose  it 
exists  and  is  located  at   z(t).   Then  it  follows  as  above 

44 


that 

lim   d(x+z(t),t)  =  s„(x), 

|t|-^GO  2 

uniformly  on  bounded   x-intervals. 

A  careful  examination  of  equation  (2.28)  satisfied 
by   d   shows  that   d   does  indeed  have  another  local 
maximum  for   |t|   large;   this  proves  the  theorem  in 
case  b) . 

Case  c)  can  be  analyzed  in  a  similar  fashion. 
The  main  difference  is  that  the  function  R(m)   vanishes 

at   m  =  m-,  -m^;      it  is  easy  to  show  that  this  causes 

1/2 
solutions  of  m,  =  +  R  '^  (m)   in  the  range 

m, -m„  <_  m(t)  <  m   to  behave  in  the  following  fashion: 
There  is  a  value   t    such  that  iT^lt^)  =  m-j^-m^. 


and  m  satisfies 


-  R-^/^(m)   for   t  <  t^ 


m,  = 
^      ,1/2 


R^/  (m)   for   t^  <  t 


o 


In  particular  m(t)   tends  to  m,   as   |t|   tends  to   oo  . 
As  in  case  b)  one  can  show  that  for   |t|   large 

there  is  exactly  one  absolute  and  one  local  maximum,  and 

that  the  asymptotic  relation  (2.31)  holds.   This  completes 

the  proof  of  theorem  2.1. 

The  time  history  of  the  local  maxima  is  as  follows: 
As   t   goes  from   -co   to   oo ,  the  amplitude  of 
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the  smaller  solitary  wave  increases  until  it  reaches  the 
value  n   ,      at  which  point  the  local  maximum  disappears. 
In  the  meanwhile  the  amplitude  of  the  larger  solitary 
wave  decreases  steadily;  in  case  c)  this  amplitude  reaches 
its  minimum  value   m-, -m^  at  some  time   T  and  some  point   X. 
The  double  wave  is  symmetric  with  respect  to  this  occurrence, 
i.e. 

(2.38)  d(X-x,T-t)  =  d(x,t)  . 

In  case  b)  the  amplitude  of  the  larger  solitary  wave 
decreases  until  it  reaches  the  value   n, ,   at  which  point 
the  local  maximiim  disappears.   Before  this  happens  however 
another  local  maximum  is  created  which  starts  increasing. 
Denote  by  T   the  time  when  the  two  maxima  are  equal,  and 
denote  by  X  the  midpoint  between  the  two  local  maxima; 
the  double  wave  satisfies  the  symmetry  relation  (2.38). 

Speaking  qualitatively  we  might  say  that  in  cases 
b)  and  c)  the  big  wave  first  absorbs,  then  reemits  the 
small  wave,  and  that  in  case  b)  the  absorption  of  the  small 
wave  raises  a  secondary  peak  on  the  big  wave. 

Up  to  a  certain  point  one  can  give  a  similar  analysis 
of  N-tuple  waves.   That  is, solutions  of  the  KdV  equation 
which  as   t  — >  +  00   split  apart  into  a  superposition 
of  n-solitary  waves  with  speeds   c,  ,  Cp,  .  .  .  c,,.   Using  the 
first   N   of  the  sequence  of  integrals  constructed  in  [5] 
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one  can  derive  an  ordinary  differential  equation  of  order 
2N  with  respect  to   x.   An  elegant  argument  of  Gardner 
(personal  communication)  shows  that  the  generalized   KdV 
operators   K  ,   n  =  1,...,  N  are  integrating  factors 
for  this  differential  operator.   In  this  fashion  we  can 
obtain  by  elimination  differential  equation  of  order   N 
for  the   N-tuple  wave,  but  the  resulting  equation  is  too 
complicated  to  yield  any  useful  information. 
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